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The neutron capture cross section of 26 Mg was measured relative to the known gold cross section 
at thermonuclear energies using the fast cyclic activation technique. The experiment was performed 
at the 3.75 MV Van-de-Graaff accelerator, Forschungszentrum Karlsruhe. The experimental capture 
cross section is the sum of resonant and direct contributions. For the resonance at £? n ,iab = 220 keV 
qq , our new results are in disagreement with the data from Weigmann et al. An improved Maxwellian 

■ averaged capture cross section is derived from the new experimental data taking into account s- and 

^\ I p-wave capture and resonant contributions. The properties of so-called potential resonances which 

influence the p-wave neutron capture of 26 Mg are discussed in detail. 
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PACS numbers: 25.40.Lw, 24.50.+g 



I. INTRODUCTION 



Neutron capture processes of neutron-rich light isotopes play an important role in astrophysical scenarios ranging 
from the so-called inhomogeneous Big Bang models to nucleosynthesis in stellar helium and carbon burning stages. 
In Inhomogeneous Big Bang Models the high neutron flux induced primordial nucleosynthesis bridges the mass 5 
and mass 8 gap |lj . Subsequent neutron capture processes of neutron-rich isotopes may even lead to a primordial 
r-process ||||. The efficiency for the production of heavy elements in such a scenario depends sensitively on the 
£^ ■ respective neutron capture rates for these light isotopes. Therefore the neutron capture cross sections have to be 
OO \ determined over a wide energy range up to about 1 MeV. 

. In massive Red Giant stars magnesium isotopes are mainly the products of hydrostatic carbon and neon burning. 

25 Mg and 26 Mg have also appreciable abundances in the outer carbon layer as a result of the reactions 22 Ne(a,7) 26 Mg, 
, 22 Ne(a,n) 25 Mg, and 25 Mg(n,7) 26 Mg jjj. He-shell burning in low mass asymptotic giant branch (AGB) stars has been 
1 1 ■ proposed as the site for the main component of the s-process ||^|. In AGB stars 26 Mg is likely to be made by 
22 Ne(cv,7) 26 Mg, as well as by 25 Mg(n,7) 26 Mg, by the decay of 26 Al, and by 26 Al(n,p) 26 Mg with radioactive 26 Al being 
synthesized through the previous H-burning shell. For the above discussed stellar scenarios neutron capture rates 
need to be known in the energy range between about 5 and 200 keV. 

The reaction rate for 26 Mg(n,7) 27 Mg that is being investigated in this work is low through its small cross section. The 
reaction 26 Mg(n,7) 27 Mg is of interest in stellar nucleosynthesis, because (i) 26 Mg is one of the most abundant isotopes 
in the cosmos, (ii) all Mg-isotopes show a strong deviation from the normal 1/v-behavior of their cross sections, and 
(iii) of the interplay between neutron production and captures for Mg-isotopes, and with the concurrential radioactive 
26 A1. 

The cross section for neutron capture processes is dominated by the nonresonant direct capture process (DC) and 
by contributions from single resonances which correspond to neutron unbound states in the compound nucleus (CN) . 
For calculating the different reaction contributions we used a simple hybrid model: the nonresonant contributions 
were determined by using a direct capture model, the resonant contributions were based on determining the resonant 
Breit Wigner cross section. In the case of broad resonances additional interference terms have to be taken into account 
which were neglected in the present work because the relevant resonances are relatively narrow. 

In this work we also discuss the effects of potential resonances in direct capture. Potential resonances are well known 
in nuclear processes (see Ref. |?J and references therein), but their influence in numerical direct-capture calculations 
has to our knowledge not been investigated before in great detail. The reaction 26 Mg(n,7) 27 Mg provides an excellent 
example for discussing the artifacts such resonances can produce in direct-capture calculations. 

Up to now the neutron capture cross section of 26 Mg was measured at thermal energies in several experiments 
(see Refs. pHLCJ], and references therein). However, at thermonuclear energies only the resonance properties of four 
resonances in^Mg(n,7) 27 Mg were determined by Weigmann et al. |[TT[] , and in that work the Maxwellian averaged 
capture cross section <a> = <a-v> / vt, with vt = ^j2kgT / /i, was calculated by the sum of the s-wave 
contribution (extrapolated from the thermal capture cross section assuming the usual 1/v law) and of the contributions 
from the four measured resonances. The p-wave contribution of the DC was neglected in that work. 
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In the following we will present our experimental procedure and results in Sect. ||. In Sect. Ill we analyze our new 
experimental data using the DC model and we discuss the properties of so-called potential resonances, in Sect. |v| we 
show the improved Maxwellian averaged capture cross section < a >, and the results are discussed in Sect. |v|. 



II. EXPERIMENTAL SETUP AND RESULTS 



The experiment was performed at the 3.75 MV Van-de-Graaff accelerator of the Forschungszentrum Karlsruhe. 
The residual nucleus 27 Mg decays with a half-life of T 1/2 = 9.458 min to 27 A1. Two 7 rays with E 1 = 843.76 keV 
(branching 71.8 %) and E 1 = 1014.44 keV (branching 28.0 %) can be detected Q. The decay properties of 27 Mg 
and 198 Au (from the neutron capture of 197 Au) are listed in Table ffl. 



A. Fast cyclic activation technique 



For residual nuclei with short half-lives the fast cyclic activation technique was developed at the Forschungszentrum 
Karlsruhe |13|. In the 26 Mg experiment the sample is irradiated for a period t^, after this irradiation time the sample 



is moved to the counting position in front of a high-purity germanium (HPGe) detector (waiting time i Wl ). The 
7-rays following the /3-decay of 27 Mg are detected for a time interval t c , and finally the sample is moved back to the 
irradiation position (waiting time t W2 ) . The whole cycle with duration T = % + t Wl +t c + t W2 is repeated many times 
to gain statistics. The 26 Mg(n,7) 27 Mg experiment was performed partly with ib = 49.6 s, t Wl = 0.4 s, t c = 49.6 s, 
i W2 = 0.4 s, T = 100 s, partly with t h = 119.4 s, t Wl = 0.6 s, t c = 119.4 s, t W2 = 0.6 s, T = 240 s. 

The accumulated number of counts from a total of N cycles, C = Y27=i where Cj, the counts of the i-th cycle, 
are calculated for a chosen irradiation time, t^, which is short enough compared with the fluctuations of the neutron 

flux, is in 



1 ( \ + \ 

C = e 7 A' 7 / 7 A" 1 [l - exp(-Atc)] exp(-At Wl ) - - ° XP , Na y [l - f h exp(-AT)] ^ 



1 - cxp(-AT) 

with 



i=i 



(2.1) 



Sr = i^exp[-(n-QAT] 

/b = EL^ • (2 - 2) 

The following additional quantities have been defined: e 7 : detection efficiency, if 7 : 7-ray absorption, / 7 : 7-ray 
intensity per decay, A: decay constant, N: the thickness (atoms per barn) of target nuclei, er 7 : the capture cross 
section, the neutron flux in the i-th cycle. The quantity /b is calculated from the registered flux history of 
a Li glass monitor which is mounted at a distance of 91 cm from the neutron production target. A simple n-7 
discrimination was applied to the slow output of the 6 Li monitor: the amplitude of the neutron signal from the 
6 Li(n,a) 3 H reaction is much larger than the amplitude from a typical 7 ray event that cannot deposit its full energy 
in the 1 mm thin 6 Li glass scintillator. 

The experimental uncertainties are reduced to a great amount by measuring the 26 Mg capture cross section relative 
to the well-known gold standard Jl3]. For this purpose our magnesium samples were sandwiched between two thin 
gold foils. Then the capture cross section of the sample er 7 is given relative to the capture cross section of the gold 
reference by: 

s _ R C s N K e*K*f*\ s ll - exp(-\ K t c )} exp(-A R f Wl )[l - cxp(-A R t b )][l - cxp(-A s r)][l - / b R cxp(-A R T)] 
°i ~ a ~t c R iVS e SKS/ 7 S A R[l - exp(-ASi c )] exp(-AS* wl )[l - exp(-ASt b )][l - exp(-A R T)][l - exp(-AST)] 

(2.3) 



B. Samples 

Our samples consisted of MgO powder enriched in 26 Mg by 98.79 % and of metallic magnesium of natural isotopic 
composition. The MgO powder was pressed to two small tablets with a diameter of d = 6 mm and a thickness of 
about 0.5 mm and 1.5 mm, respectively. The MgO tablets were put into a thin Al foil to ensure that no powder is lost 



2 



during the measurements. The mass of the tablets had to be determined very carefully because the MgO powder is 
very hygroscopic. Therefore, the tablets were heated at T w 1000° C for several hours, and then they were weighted 
during the cooling phase. We removed the tablets from the oven at about 300° C. After this procedure we found that 
the measured weight of the tablets was stable for a few minutes, but then the weight began to increase by about 10 
to 20 percent within a few hours. 

As a further check we repeated one activation measurement with a sample of metallic nat Mg with the same geometric 
shape (diameter 6 mm, thickness 1 mm, mass 47.77 mg). The result agreed perfectly with the cross section measured 
with the enriched samples. This confirms that the weighting procedure worked well, that the enrichment of 26 Mg in 
the MgO samples was correct, and that the influence of degraded neutrons scattered by hydrogen in the hygroscopic 
MgO samples is negligible. The properties of the samples are summarized in Table ||. 



C. Neutron production, Time-of-Flight measurements 

The neutrons for the activation were generated using the 7 Li(p,n) 7 Be reaction. We measured the capture cross 
section in the astrophysically relevant energy region of several keV at six different neutron energies which are obtained 
by different combinations of proton energy and thickness of the lithium target: (?) Close above the reaction threshold 
at B Pi iab = 1881 keV one obtains a quasi-Maxwellian spectrum with ksT = 25 keV |15j if the energy loss of the 
protons in the lithium layer of the target reduces the proton energy below the reaction threshold, (ii) Very close to 
the reaction threshold a narrow spectrum with 12 keV < -E n ,iab < 52 keV was obtained. Note that in these cases (i) 
and (ii) the neutron energy spectrum is independent of the thickness of the lithium layer. We used a lithium layer 
of 30 fim which is by far thick enough for that purpose. (Hi) Using four different higher proton energies we obtained 
four different neutron energy spectra with higher neutron energies and an energy width of about 20-50 keV. Now the 
energy width depends sensitively on the thickness of the lithium layer which was about 2.6 /zm in our experiment. 
The thicknesses of the different thin lithium targets were determined by two independent methods: first during the 
evaporation of the metallic lithium layer and second from the width of the time-of- flight spectra (see below) . Both 
methods give the same results within an uncertainty of about 10%. 

For the activation runs beam currents in the order of 80 to 100 /xA were used. Because of the electric power of the 
proton beam of about P«2 MV • 100 /iA s» 200 W the copper backing of the thin lithium layer was water-cooled, 
and the temperature of the copper backing remained below 85°C during the whole experiment. 

The neutron energy spectra were controlled by time-of-flight (TOF) measurements using the pulsed beam which 
is available from the Van-de-Graaff accelerator. Typically the repetition frequency is 1 MHz, and the pulse width is 
about 10 ns. The beam current in the pulsed mode was about 1 to 4 /iA. The neutrons were detected in the 6 Li glass 
monitor, and the TOF of the neutrons was measured between the fast timing output from the 6 Li monitor (start) 
and the delayed pickup signal from the pulsed beam (stop). Especially at energies close to the resonances in the 
26 Mg(n,7) 27 Mg reaction the energy of the neutrons has to be determined very accurately. A typical TOF spectrum 
at i?iab = 0° is shown in Fig. [jj 

The result of our experiment is an average capture cross section 

a = J a(E nMh )$(E nMh )dE nMh (2.4) 

where the neutron energy distribution ^(-En.iab) is normalized to 1: 

J ®(E nMh )dE nAah = 1 (2.5) 

The neutron energy distribution $(£ , n ,iab) was calculated from the measured TOF distribution <i>(T) by the relation 
&(T)dT = $(£ , n ,iab)rf^' ni iab- Additionally several facts have to be taken into account: (?) the angular dependence 
of the 7 Li(p,n) 7 Be reaction cross section which was measured at several angles using the 6 Li monitor |l6| , (ii) the 
geometry of the neutron production target and the sample, (Hi) the energy loss of the protons in the lithium layer with 
a thickness d = 2600 nm± 10% |17 , (iv) the low-energy tail of the neutron energies in the TOF spectrum (see Fig. [I]). 



(Hi) and (iv) are relevant only for the measurements at higher neutron energies. The neutron energy distribution 
with -En, lab = 208.3 keV is shown in Fig. |^. It has to be pointed out that it is difficult to measure the effective neutron 
energy distribution $(£'„, lab) directly because of the extended lithium target and the extended sample. However, all 
ingredients (i)— (iv) for the calculation of $(£ , n! i a b) were determined experimentally, and so the uncertainties of the 
calculation are small. 
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D. 7 ray detection 



The 7 rays from the decay of 27 Mg and 198 Au were detected using a high purity germanium (HPGe) detector with a 
relative efficiency of 60 % (compared to a 3" x 3" Nal detector). The energy-dependent efficiency of the HPGe detector 
was calculated using the detector simulation program GEANT pl| . The simulation calculation was controlled by an 
experimental determination of the efficiency with several radioactive calibration sources. The results for the relative 
efficiencies are 

e f 1L80keV ; = 1.752 ± 2% and b^J^D. = 2 . 016 ± 2% for d = 30 mm 
e 7 (843.76 keV) e 7 (1014.44 keV) 

^{ 41L8 ° keV ; = 1.746 ± 2% and hj^Ml = 2 . 007 ± 2% for d= 35 mm 
e 7 (843.76 keV) e 7 (1014.44 keV) 

for the close geometry used in our experiment (the distance between sample and detector was 30 mm for the fcgT = 
25 keV runs and for the run very close above the threshold, and 35 mm for the other runs). Additionally, the 7-ray 
absorption in the sample and in the gold foils was calculated using the cross section tables provided by National 
Nuclear Data Center, Brookhaven National Laboratory, via WWW, and based on Ref. j[9). This small correction is 
in the order of a few %. A typical spectrum from the HPGe detector is shown in Fig. |^. 



E. Experimental Results 



The new experimental data are listed in Table [II, and they are compared to a DC calculation (see Sect. [II) in 
Fig. ^. Two of the four known resonances at E nMh = 68.7 keV, 220 keV, 427 keV, and 432 keV |llj] have influence 
on our measured cross sections at fc^T = 25 keV (resonance at 68.7 keV) and at -E n .iab = 178 keV and 208 keV 
(resonance at 220 keV). The two resonances above 400 keV do not influence our measurement because they are 
narrow (r < 2 keV) |J. 

The influence of a narrow resonance at the energy E rcs on the measured capture cross section can be calculated by 

^rcs = / <WE c .m.)*(£c.m.)dEcm. ~ ^nTab) -^±i_ 2 7T — (yry) (2.6) 

where asw{E c ,m.) is the Breit-Wigner shaped capture cross section in the resonance, $(E) the neutron energy 
distribution, the wave number at the resonance energy E Ies = h 2 {kl c ^) 2 /2n in the center-of- mass system, and 
luj is the resonance strength. 

The resonance strength of the -Eii,iab = 220 keV resonance was measured by Weigmann et al. [|llf| : = 3.75 ± 
0.16 eV. This resonance gave the dominating contribution to the capture cross section for the neutron energy spectrum 



^K-Enjab) shown in Fig. g (corresponding to E® lab = 2003 keV). ^From Eq. 2.i one obtains a les — 1.78 ±0.27 mb. In 
contradiction to that result our measured capture cross section is only a — 0.735 ± 0.040 mb as shown in Fig. |^! 

Our method for measuring neutron capture cross sections is mainly tailored for non-resonant capture. Nevertheless, 
our experimental result definitely excludes the result of Weigmann et al. jllj for the resonance at 220 keV. To estimate 
the uncertainties of our calculated <f>(-E n jab) we calculated energy distributions $(-E n ,iab) also at proton energies 
-E"p,iab = £p iab i 2 keV which is much larger than the uncertainty of the accelerator calibration which was confirmed 
by the TOF measurements. These results for £ P) i a b = E® lab ± 2 keV are also shown in Fig. |2| as dashed lines. The 
resulting uncertainties for $(£ , i a b = 220 keV) and a res are about 15%. 



If one subtracts the expected non-resonant DC contribution cjdc = 0.095 ± 0.010 mb (see also Sects. Ill A and IV) 
from our measured capture cross section at £ n ,iab = 208.3 keV (-E c .m. = 200.3 keV) we obtain a resonant contribution 
of <t = 0.64 ± 0.05 mb, which is a factor of about 2.8 lower than the result from Weigmann et al. JLT| |. Therefore, we 
obtain a resonance strength of uj^ = 1.34 ± 0.24 eV for the -E n ,lab = 220 keV resonance. Interference effects between 
the DC and the -En, lab = 220 keV resonance can be neglected because the total width of this resonance is relatively 
small: T = 0.44 ± o!o7 keV (Xl| - 

A similar analysis of the 68.7 keV resonance has relatively large uncertainties because only the tail of the neutron 
energy spectrum with ksT — 25 keV overlaps with this relatively weak resonance (ujj = 73 ± 20 meV from Ref. [ pT| , 
at least more than one order of magnitude smaller than the -E n .i a b = 220 keV resonance), and the dominating 
contributions to this data point at fcgT = 25 keV come from s- and p-wave capture. 
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III. THEORETICAL ANALYSIS 



A. Direct capture model and folding potentials 



The theoretical analysis was performed within the Direct Capture (DC) model. The DC cross section is given by 
©0 



a 1 % 



Ma M a Mb , 7r-v 



The quantities I a, Ib and S a (Ma, Mb and M a ) are the spins (magnetic quantum numbers) of the target nucleus 
A, residual nucleus B and projectile a, respectively. The reduced mass in the entrance channel is given by /i. The 
polarization 7r 7 of the electromagnetic radiation can be ±1. The wave number in the entrance channel and for the 
emitted radiation is given by k a and fe 7 , respectively. The transition matrices T = T E1 + T E2 + T M1 depend on the 
overlap integrals 

%££ M . = J ^u hjblB (r)0^(r)xi aja (r) (3.2) 

The radial part of the bound state wave function in the exit channel and the scattering wave function in the entrance 
channel is given by Ui b j b i B (r) and xi a j a ( r )-> respectively. The radial parts of the electromagnetic multipole operators 
are well-known. The calculation of the DC cross sections has been performed using the code TEDCA p2|| . 

For the calculation of the scattering wave function x(r) and the bound state wave function u(r) we used systematic 
folding potentials: 

V{R) =\V F (R) =X J y > p a (r 1 V A (r 2 )« cff ( J E,p a ,p yl , s )dr 1 dr 2 (3.3) 

with A being a potential strength parameter close to unity, and s = |R+r2 — ri | , where R is the separation of the centers 
of mass of the projectile and the target nucleus. The density of 26 Mg has been derived from the measured charge 
distribution P3fl , and the effective nucleon-nucleon interaction v e s has been taken in the DDM3Y parametrization 
p4j . The resulting folding potential has a volume integral per interacting nucleon pair Jr = 500.41 MeVfm 3 (A = 1) 
and an rms-radius r rms = 3.676 fm. Details about the folding procedure can be found for instance in [ps| , the folding 
potential has been calculated by using the code DFOLD [^6| . The imaginary part of the potential is very small be cause 
of the small flux into reaction channels and can be neglected in our case. This fact will be discussed in Sect. [II D| . 



The potential strength parameter A can be adjusted either to the binding energy Eb (bound state wave function) or 
to the scattering length (s-wave scattering at thermal energies). 

The Pauli principle was taken into account for the bound state wave functions by the condition 

Q = 2N + L (3.4) 

where Q, N, and L are the number of oscillator quanta, nodes in the wave function, and angular momentum number. 
For even-parity states one has to take Q = 2 in the sd shell, and for odd-parity states Q = 3. The parameters A for 



several bound state wave functions are listed in Table [V. Note that the A values do not differ more than 30% from 



unity for all bound states (A ~ 0.9 for Q — 2 states with positive parity, A 1.25 for Q — 3 states with negative 
parity). 

The total capture cross section is given by the sum over all DC cross sections for each final state i, multiplied by 
the spectroscopic factor C 2 S which is a measure for the probability of finding 27 Mg in a ( 26 Mg ® n) single-particle 
configuration 

<7 7 =^C 2 S iCT ° c . (3.5) 

i 

For the 26 Mg(n,7) 27 Mg reaction one has to take into account all final states below the neutron threshold of 27 Mg. 
However, in practice only El transitions to final states with large spectroscopic factors are relevant in the keV energy 
region. Additionally, at thermal energies Ml and E2 transitions contribute to about 17% to the thermal capture cross 
section 
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B. Spectroscopic factors 



For the theoretical prediction of the neutron capture cross section the knowledge of the spectroscopic factors 
(SF) is essential. These SF have been determined from the DWBA analysis of different transfer reactions like 
e.g. 26 Mg(d,p) 27 Mg. Unfortunately, the experimental SF from the different experiments only agree within a fac- 
tor of 2 (although these experiments claim uncertainties in the order of 10-20%). Additionally, SF from theoretical 
shell model calculations are in average roughly a factor of 2 smaller than the experimental SF. These discrepancies 
are listed in Table IV . 

To reduce these uncertainties the following procedure was applied. At thermal energies SF can be determined with 
high accuracy from the ratio of the experimental capture cross section to the calculated DC cross section. The DC 
cross section can be calculated with small uncertainties because the uncertainties in the calculation of the scattering 
wave function can be minimized by adjusting the potential strength (parameter A) to the neutron scattering length 
@. In the case of 26 Mg the scattering lengths & b ound = 4.89 ± 0.15 fm @ and 6 frco = (b hound - Z b nc ) A/ (A + 
1) = 4.725 ± 0.15 fm lead to A = 0.9694 ± 0.024. Using this value for A we obtain for the two dominating El 
transitions at thermal energies: C 2 S(3/2~ ,3560 keV) = 0.360 ± 0.054 and C 2 5(l/2~, 4827 keV) = 0.260 ± 0.0390. 
The calculation includes the uncertainties of the potential strength parameter A, of the thermal capture cross section 
(^thermal = 39.0 ± 0.8 mb [H), and of the branching of the thermal capture H|J. 

The result for the SF is in average 14.1% smaller than the result of Meurders et al. [^7j obtained from (d,p) 
measurements for the two odd-parity levels for which the SF can be determined accurately as described above. 
Therefore, we renormalized the experimental SF from Ref. |27| by the factor 0.852 ± 0.091. The experimental data 
of Ref. fl27| were used because only in that experiment all relevant levels were analyzed. The result for the DC cross 
section is compared to the experimental results in Fig. |[ One can clearly see the transition from the 1/u-behavior (s- 
wave capture) to the u-behavior (p-wave capture) at E « 20 keV. This transition is typical especially for neutron-rich 
nuclei in the sd shell with many bound states with even parity and large spectroscopic factors. 

The overall agreement between the experimental data points which are not affected by the known resonances and 
the DC calculation is quite good. However, our DC calculation overestimates the p-wave capture. There are two 
possible explanations: (i) the spectroscopic factors of the final states with even parity (which cannot be determined 
by the procedure described above) are too large, (ii) the optical potential for the p-wave differs slightly from the 
s-wave potential. 



C. Potential resonances 



In the following our discussion is restricted to folding potentials with one free parameter A. The same arguments 
also hold for Saxon- Woods or other potentials, but in many cases the discussion may be more complicated because 
the larger number of parameters of the Saxon- Woods potential. 

The real part of the optical potential for the incoming p-wave was assumed to be identical to the s-wave potential 
which was adjusted to the neutron scattering length at thermal energies. However, a weak dependence (in the order 
of about 10-20%) of the optical potential on the parity and/or the angular momentum of the partial wave was found 
in many systems (see e.g. Ref. p7|). This weak energy dependence has important consequences for the calculation 
of the neutron capture cross section of 26 Mg, because the 26 Mg-n potential shows a so-called potential resonance for 
the p-wave close to the astrophysically relevant energy region. This fact leads to a very sensitive dependence of the 
capture cross section on the potential strength (i.e. the strength parameter A). In Fig. pj we show the dependence of 
the capture cross section on the potential strength parameter A. Changing A from 0.8 to about 1.0 does not change 
the qualitative shape of the cross section curve, but the absolute value changes by a factor of 4.5 at 1 MeV. For 
A > 1.0 the resonant behavior becomes more and more visible, and for A = 1.10 the resonance disappears below the 
threshold and leads to a bound state at Eb = —135.3 keV. 

Such potential resonances are characterized by a large width because the simple two-body potential model assumes 
full single-particle strength for such resonances. For very broad resonances (r E) it is difficult to give precise 
numbers for the position and the width of a resonance (see e.g. the discussion in Ref. j38[). To avoid such problems 
of definition we show the p-wave phase shift Si = i(E) for several parameters A in Fig.Tfn and we define the energy 



1 The numerical values are slightly different from a previous work |35| because of the improved experimental data |8j compared 
to Q and because of an improved adjustment of the potential strength to the free coherent scattering length 6f ree instead of 
the bound frbound- 
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of the resonance as the maximum slope of the phase shift, and the width is given by the slope at that energy: 
r = 2 [(dSi = i(E) / dE)\E=E IOS ]~ 1 ■ The resulting values for E ICS1 V, and the phase shift 5i=i(E Ies ) are listed in Table [V|. 

For a proper adjustment of the p-wave potential strength the knowledge of the experimental scattering phase 
shifts is desirable. However, to our knowledge these phase shifts have not been measured. Nevertheless, there is an 
experimental hint that such a potential resonance is not only an artifact of the theory, but exists in nature: a broad 
structure was observed in a measurement of the total 26 Mg-n cross section at about E 300 keV At that energy 
our potential model predicts a width of several hundred keV. 

In general, several potential resonances have been proved experimentally. A typical example is the first 3 + state 
in 6 Li which can be found as a resonance in the 2 H(a,7) 6 Li capture reaction |39[ ). Because of the almost pure single- 
particle structure of 6 Li = 2 H <g> a the properties of this resonance can be predicted by the potential model. For 
heavier nuclei an almost pure single-particle structure cannot be expected but there still exist broad resonances with a 
relatively strong single-particle component. This leads to an overestimation of the width because the simple two-body 
potential model assumes full single-particle structure for potential resonances. As a consequence, the potential model 
is not able any more to predict the properties of such resonances with high accuracy; only if the spectroscopic factor 
of such a resonance is known, the width can be calculated using the relation 

T cxp = C 2 ST s . p . (3.6) 

where r s p . is the width calculated from the potential model. Furthermore, the strength of the theoretical single- 
particle resonance is fragmented into several states which can be identified experimentally. This can already be seen 
from the fact that the two bound p-wave states which dominate the thermal capture cross section of 26 Mg(n,7) 27 Mg 
have SF in the order of about 0.2 to 0.4. 

It has to be pointed out that the so-called "non-resonant" energy region can be calculated by the simple potential 
model only if the potential is carefully adjusted because the potential model always contains potential resonances, 
and as one can see from Figs. || and |], there is no clear definition for "resonant" or "non-resonant". A wrong 
adjustment may lead to (i) an overestimation of the capture cross section if the potential contains a single-particle 
resonance which does not exist in the experiment, and (ii) an underestimation of the capture cross section if the 
adjustment neglects strong (and broad) single-particle resonances. 

The success of many previous calculations using the DC model, which seems to be in contradiction to the above 
statements, can be explained by two facts: (i) in several cases the influence of potential resonances is small in the 
energy region which was analyzed in the calculation (this is especially the case when the main contribution for the 
DC transition comes from regions far outside the nucleus, see e.g. Refs. [|o| or |4]J| ), (ii) even wrong DC calculations 
can be cor rected using adjustable spectroscopic factors (uncertainties in the order of a factor of 2 are realistic, see 



Sect. [Ill B| !). 

^From this point of view the rough agreement of the calculated DC cross section with the experimental results is 
already somewhat surprising because the range of the capture cross section (see Fig. ||) covers five orders of magnitude 
when the potential strength changes by only ±15% from about 0.95 (determined from the s-wave scattering) down 
to 0.8 and up to 1.1. Such a dramatic behavior of the so-called "non-resonant" direct capture can appear in any 
DC calculation. To the best of our knowledge, these facts were not completely discussed in all the previous papers 
dealing with DC. In conclusion, we point out that the calculation of the DC has to be done with great care for each 
nucleus because otherwise deviations between the predicted and the experimental capture cross section not only by 
a factor of two, but even by orders of magnitude can appear, and both overestimations and underestimations of the 
experimental cross sections are possible! 



D. The imaginary part of the optical potential 

As pointed out by Lane and Mughabghab G3] , the capture amplitude has to be derived from the full optical 



potential, i.e. real and imaginary part (see also Eq. 3.2). However, in the case of 26 Mg(n,7) 27 Mg the imaginary part 
of the optical potential can be neglected for the following reasons. The only open inelastic channel is neutron capture. 
Then the capture cross section a 1 is directly related to the reflexion coefficient r)i by: 

V^c.m.J / 

Using the real folding potential with A = 0.9694 (from Sect, [ill B|) , an imaginary Saxon- Woods potential with depth 

1/3 

Wo, reduced radius Rq = 1.3 fm (R — Ro ■ Aj, ), diffuscness a = 0.65 fm, and the experimentally determined values 
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of 07, one obtains a depth of the imaginary part in the order of Wq « 100 eV for the optical s- and p-wave potential. 
This very weak imaginary part does not significantly influence the calculated DC cross section. 

A significant reduction of the DC cross section is only obtained when Wq is increased by more than three orders of 
magnitude. This means that the calculations of the DC cross section in Fig. || are practically not influenced by the 
imaginary part of the potential. Only potential resonances with very huge capture cross sections (with A ~ 1.08) will 
be damped slightly by the imaginary part of the potential. 



IV. ASTROPHYSICAL REACTION RATES 



The astrophysical reaction rate Na < u ■ v > can easily be derived from the Maxwellian averaged capture cross 
section <(i>~<a-v> / vt- The Maxwellian averaged cross section of the reaction 26 Mg(n,7) 27 Mg was calculated 
from the following three ingredients. First, the thermal s-wave capture taken from Ref. |3q | was extrapolated to the 
thermonuclear energy region using the well-known 1/v law. This 1/v behavior is reproduced by the DC calculations. 
Therefore, the uncertainty of the s-wave contribution is below 5%. Second, the p-wave capture cross section which 
dominates at thermonuclear energies is roughly proportional to the velocity v: 

a = C-VE. (4.1) 



The constant C was determined from the experimental data: C = 180±30 fib/y/MeV. Because only three experimental 
cross section data are not influenced by resonances, and because one has to subtract the s-wave contribution from the 
experimental capture cross section, we estimate an uncertainty of less than 20% for the p-wave contribution. Note 



that the the p-wave contribution of the DC for the data point at E n .i a b = 208.3 keV was also taken from Eq. 4.1, 
and the s- and p-w ave c ontribution were subtracted to determine the resonance strength of the -En, lab = 220 keV 
resonance (see Sect. |lIEj). Third, the resonance contributions were reduced by a factor of 2.8 compared to Ref. Jll]] 
(see discussion in Sect. [II E| ). This factor was only determined for the En, lab = 220 keV resonance, and the uncertainty 
of the sum of the resonant contributions is roughly a factor of 2. 

Finally, this leads to a Maxwellian averaged capture cross section < a > as shown in Fig. fj]. At low energies 
(ksT < 10 keV) s-wave capture is dominating, and from about ksT « 20 keV the p-wave capture exceeds the 
s-wave contribution. The resonances become dominant at higher energies (ksT « 50 keV). By accident, our result 
is relatively close to that of Ref. O): the sum of our measured p-wave capture and the by a factor of 2.8 reduced 



resonant contributions is somewhat larger at ksT < 30 keV than the calculation of Ref. |11| which neglected the 
p-wave contribution, and at higher energies (ksT > 40 keV) the result of Ref. [[flj is higher than our new result due 
to the overestimation of the resonance strengths in that work. 



V. SUMMARY AND CONCLUSION 



There are three main results of this work: (i) the DC cross section of the reaction 26 Mg(n,7) 27 Mg was measured 
for the first time between the known resonances at £" n ,iab = 68.7 keV and E n .i a b = 220 keV , (ii) the properties of 
potential resonances are discussed in detail for the p-wave capture of 26 Mg(n,7) 27 Mg, and (Hi) the resonance strength 
of the i? n , lab = 220 keV resonance was determined to be uj"f = 1.34 ± 0.24 eV by a careful analysis of the neutron 
energy spectrum $(E n ,i a b)- This last result is in clear contradiction to the results of Ref. |tl|. Unfortunately, in 
Ref. pHl the experimental information on the capture measurements is very limited, but it is stated in that work that 
the enriched 26 Mg- and 25 Mg-samples were analyzed in the same way. This fact makes the experimental results of 
Ref. jll| for the capture reaction 25 Mg(n,7) 26 Mg at least questionable. A strongly reduced 25 Mg(n,7) 26 Mg reaction 
rate might have influence on the s-process nucleosynthesis because of the role of 25 Mg as neutron poison. In this 
paper we do not want to speculate about astrophysical consequences of such a reduced reaction rate; instead, we think 
that primarily a new measurement of the resonance properties of 25 Mg is necessary to obtain reliable reaction rates. 
Unfortunately, the reaction 25 Mg(n,7) 26 Mg cannot be analyzed using our activation technique because the residual 
nucleus 26 Mg is stable. 
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TABLE I. Decay properties of the residual nuclei Mg and Au. 



Residual 


T\/2 


E-, 


Intensity per decay 


nucleus 




(keV) 


(%) 


"'Mg 


9.458 ± 0.012 min 


843.76 


71.8 ± 0.4 






1014.44 


28.0 ± 0.4 


198 Au 


2.69517 ± 0.00002 d 


411.80 


95.50±0.096 



TABLE II. Properties of the magnesium samples and the gold foils. 



Isotope 


Chemical 


Isotopic 


mass 




form 


composition (%) 


(mg) 


2b Mg 


MgO 


98.79 ± 1.90 (enriched) 


29.20 ± 0.05 


26 Mg 


MgO 


98.79 ± 1.90 (enriched) 


99.50 ± 0.20 


26 Mg 


metallic 


11.01 ±0.02 (natural) 


47.77 ± 0.02 


197 Au 


metallic 


100 (natural) 


15.8 - 16.5 



TABLE III. Experimental capture cross section (j CX p of the reaction 26 Mg(n,7) 27 Mg compared to the calculated DC cross 
section cr 7 calculated with A = 0.9694 (see Sect. [llj| and Fig. ^|). 

-En. lab (keV) Q- C x P (n,7) (pb) a 7 (pb) 

kT = 25.3 meV (39.0 ± 0.8) mb a 32.5 mb 

kT = 25 keV 124.2 ± 7.5 92.4 

33 56.8 ±3.0 98.1 

102 82.1 ±4.5 151.2 

151 84.9 ±5.0 189.9 

178 188 ±28 210.0 

208 735 ± 40 236.4 

a from Rcf. li 



TABLE IV. Spectroscopic factors of Mg = Mg gi n taken from different experiments ]27|-|3l| and from shell-model 
calculations p^ , ^2|j33| . Note that the spectroscopic factor of the state at E x = 4149.8 keV was determined only in one 
experiment. Ad dition ally, the potential strength parameter A is given, which is adjusted to the binding energy of each bound 
state (see Sect. [[II A ). 



E x (keV) r C 2 5- cxp C'ffth C 2 S (nn) X_ 

T72 T 0.44 - 1.07 0.43 - 0.70 - 0.967 
984.7 3/2+ 0.37- 0.80 0.28 - 0.45 - 0.947 
1698.0 5/2+ 0.13 - 0.31 0.02 - 0.14 - 0.922 
3559.5 3/2" 0.34 - 0.56 - 0.360 ± 0.054 1.280 
3760.4 7/2" 0.40 - 0.80 - - 1.360 
4149.8 (5/2+) 0.03* 0.001 - 0.06 - 0.827 
4827.3 (l/2~) 0.32 - 0.39* - 0.260 ± 0.039* 1.208 

* assuming (J n ) as given in column 2 
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TABLE V. Calculated resonance energies E rcB and widths T of a p-wave potential resonance for potential strength parameters 
A from 0.90 to 1.10. For A < 0.9 the potential resonance becomes very broad. 



A 


E TCB (keV) 


T (keV) 


5i=i(E ieB ) 


0.90 


795 


4130 


16° 


0.92 


755 


3170 


20° 


0.94 


705 


2400 


24° 


0.96 


647 


1785 


28° 


0.98 


562 


1284 


32° 


1.00 


485 


878 


38° 


1.02 


381 


553 


44° 


1.04 


274 


300 


51° 


1.06 


156.5 


116 


62° 


1.08 


28.6 


8.2 


78° 


1.10 


-135.3 a 







(deg) 



a bound state with Q = 2N + L = 3, N = 1, L = l 
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800 



1000 



FIG. 1. Typical neutron TOF spectrum from a thin lithium target (thickness 2600 nm ± 10%) at ■& — 0° and at a distance 
of D = 91.3 cm. One can see the narrow peak from prompt 7 rays at channel 932. The broad neutron peak starts at channel 
702 (£™f a x b = 134.5 keV), and the tail of the neutron energy distribution ends at channel 615 (£™j° b = 71.5 keV). The bar 
over the neutron peak indicates the expected energy width calculated from the energy loss of the protons |l7j and from the 
reaction kinematics neglecting the low-energy tail from the diffusion of lithium into the copper backing. The time resolution is 
about 10 ns (FWHM of the 7 peak). 




150 160 170 180 190 200 210 220 230 240 250 
E n ,iab (keV) 



FIG. 2. Calculated neutron energy distribution <I?(-En,iab) with E n = 208.3 keV which was used for the measure- 
ment of the -E n , lab = 220 keV resonance jllj. For comparison we also show the energy distributions which are ob- 
tained from a by 2 keV increased/decreased primary proton energy (dashed lines). The value at the resonance energy is 
$(£ n ,i ab = 220 keV) = 22.8 MeV -1 ± 15% (indicated by the data point). 
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4000 
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8000 



FIG. 3. Energy spectrum of the HPGe detector measured during the activation of Mg. In the insets the relevant areas 
around 411.80 keV (decay of 198 Au), and 843.76 and 1014.44 keV (decay of 27 Mg) are shown. The main background lines come 
from the decay of 56 Mn (£ 7 = 846.76 keV, only 3 keV above the 27 Mg decay line), from 28 Al, and from 122 ' 124 Sb. Note the 
logarithmic scale of the full spectrum and the linear scale of the three insets. 
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E, m . (MeV) 

FIG. 4. Experimental capture cross section of the reaction 26 Mg(n,7) 27 Mg compared to a DC calculation. All circles are 
the result of this experiment. The full circles (•) show the non-resonant part, the open circles (o) are influenced by the known 
resonances at B n ,lab = 68.7 keV and S n ,i a b = 220 keV JTTJ which are indicated by the vertical arrows. The thermal capture 
cross section was taken from Ref. ||. The horizontal error bars show the FWHM of the neutron energy spectra. Note the 
logarithmic scale of the big figure and the linear scale of the inset; the data point at -E n ,iab = 208.3 keV does not fit in the 
linear inset. 
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I 1 I 1 
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with C S(g.s.) = 1 




0.6 

,,, (MeV) 
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FIG. 5. Calculated capture cross section of the ground state transition 26 Mg(n,7) 27 Mg g . s . with C7 2 5(g.s.) = 1.0 (El 
transition from the incoming p-wave to the bound s-wave). The capture cross section depends sensitively on the potential 
strength parameter A of the incoming p-wave which was varied from A = 0.80 (lowest capture cross section) up to A = 1.10 
in steps of 0.02 (from bottom to top, exceptions are indicated in the figure). Note that the capture cross section varies up 
to five orders of magnitude when the potential strength is increased/decreased by about 15% from the value for the s-wave 
(A = 0.9694)! 
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FIG. 6. Scattering phase shifts Si-i(E) calculated for different potential strength parameters A from 0.80 to 1.10 in steps 
of 0.02 (from bottom to top). Note that the resonance position is shifted towards lower energies as the potential strength 
increases. For A = 1.10 the resonance disappears below the threshold (see also Table [yj). 
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FIG. 7. 
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Maxwellian averaged capture cross section <<7> = <<t-?;> / vt, 



100 



for the reaction 26 Mg(n,7) 27 Mg. The total 



capture cross section is given by the sum of s-wave capture, p-wave capture, and resonant contribute 
also show the result from Ref. fill. 



For comparison we 
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